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EVALUATION 


This  study  was  performed  in  support  of  the  overall  program  of  the 
Solid  State  Applications  Section  directed  toward  developing  adequate 
reliability  screening  and  qualification  testing  sequences  for  micro- 
electronic devices  in  accordance  with  the  Reliability  Technical  Program 
Objective  No.  13,  It  successfully  met  its  objective  which  vjas  to  provide 
a means  for  predicting  the  stress  on  the  lid  seal  of  a hermetic  package 
under  various  levels  of  pressure  and  constant  acceleration.  These  levels 
are  included  in  the  screening  procedures  imposed  by  MIL-M-38510  and 
MIL-STD-883.  The  results  of  the  study  will  be  used  by  the  Air  Force  to 
establish  effective  test  levels  as  a function  of  package  size  and 
material  and  also  by  part  manufacturers  as  design  guidelines.  The 
equations  which  have  been  developed  will  first  be  verified  by  experi- 
mentation and  then  used  for  formulating  revised  screening  requirements 
in  MIL-STD-883A,  "Test  Methods  and  Procedures  for  Microelectronics," 
in  the  preparation  of  proposed  Method  5008,  "Test  Procedures  for  Hybrid 
and  Multichip  Microcircuits"  to  MIL-STD-883A  and  in  support  of  the  JLC 
Joint  Technical  Coordinating  Group  Subtask  on  Hybrid  Microcircuit 
Technology  Base. 

PETER  F.  MANNO 

Solid  State  Applications  Section 
Reliability  Branch 
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SUMMARY 


■ ! 

'-1 

An  analysis  Is  made  of  the  tensile  stresses  in  the  lid  seal  and 
the  lid  deflections  for  a rectangular  flat-pack  under  external  pressure. 
On  the  basis  of  this  analysis,  formulas  and  charts  are  presented  to 
facilitate  (a)  the  proper  design  of  the  package  so  that  it  will  retain 
its  hermetlcity  under  a given  screening  pressure  and  (b)  the  selection 
of  the  proper  pressure  to  use  in  the  hermetlcity  screening  of  an  already 
designed  package. 

Information  is  also  given  on  the  approximate  equivalence  of  external 
pressure  and  centrifuge  acceleration  in  regard  to  the  seal  stresses  and 
lid  deflections  of  a rectangular  flat-pack. 

INTRODUCTION 

External  pressure,  per  Method  1014  of  reference  1,  is  generally 
accepted  as  a means  of  screening  out  non-hermetic  or  potentially  non- 
hermetlc  microelectronic  packageo.  Knowing  that  this  is  the  screening 
device  to  be  employed,  the  designer  and  the  user  of  such  packages  are 
then  faced  with  the  following  two  Important  questions,  respectively: 

(a)  Given  the  screening  pressure  to  which  the  package  will  be 
subjected,  what  should  the  design  features  of  the  package 
be  in  order  that  a seal  of  good  quality  will  retain  its 
hermetlcity  under  that  pressure? 

(b)  Giver  a package  that  is  already  designed  and  available, 
what  screening  pressure  should  be  used  in  order  to  weed  out 
(i.e. , destroy  the  hermetlcity  or  aggravate  the  non-hermeti- 
city  of)  those  seals  at  the  low  end  of  the  quality  spectrum? 
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This  report  is  Intended  as  an  aid  in  answering  these  two  questions, 
with  particular  reference  to  flat-pack  configuration  shown  in  Figure  1. 
The  following  two  hypotheses  form  the  main  basis  of  the  present 

work: 

(a)  Under  external  pressure  the  lid  of  the  package  may  be 
regarded  as  a uniformly  loaded  rectangular  plate  with 
edges  elastically  restrained  against  rotation,  this 
restraint  being  furnished  by  the  walls  of  the  package. 

(b)  The  external  pressure  weeds  out  un:iatlsfactory  seals  by 
creating  tensile  stresses  in  the  sealing  material  that 
exceed  the  tensile  strength  of  that  material.  (In  poor 
quality  seals  (for  example,  those  containing  voids  or 
inclusions  of  foreign  matter)  this  tensile  strength  will 
presumably  be  lower  than  in  good  quality  seals.) 

SYMBOLS 

The  following  symbols  will  be  used  for  the  geometric  parameters 

and  elastic  constants.  Other  symbols  will  be  Introduced  as  needed. 

a =■  width  of  lid  (in.) 
b = length  of  lid  fin.)  (b^a) 
t = hlckness  of  lid  (in.) 
w = thickness  of  walls  (In.) 
h =»  Inside  height  of  walls  (In.) 

Young's  modulus  of  lid  material  (psi) 

=■  Young's  modulus  of  wall  material  (psi) 

V =»  Poisson's  ratio  of  lid  material 
c 

V »•  Poisson's  ratio  of  wall  material 
w 

D =•  E t^/[12(l-v  ^)]  “ plate  flexural  stiffness  of  lid 
^ ^ 3 ^2 

D •=  E w /[12(l-v  )]  ” plate  flexural  stiffness  of  walls 

w w w 
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ELASTIC  RESTRAINT  FURNISHED  BY  WALLS 


The  walls  will  be  regarded  as  wide  vertical  beans  of  height  h, 
with  the  local  rotation  0 (In  radians)  at  the  upper  end  proportional 
to  the  local  Intensity  M of  the  bending  moment  per  unit  width  exerted 
upon  the  wall  by  the  edge  of  the  lid.  This  relationship  can  be 
expressed  as 

M = k0  (1) 

where  k Is  a proportionality  constant.  Tlie  value  of  k depends  upon  the 
restraint  which  the  base  of  the  package  furnishes  to  the  lower  edge  of 
the  wall.  Two  bracketing  cases  can  be  envisioned  and  are  illustrated 
In  Figure  2.  In  one  of  these  (Fig.  2(a))  the  base  Is  so  stiff  in  bending 
(perhaps  because  of  a thick  substrate  bonded  to  its  entire  surface)  that 
it  prevents  essentially  any  rotation  of  the  lower  edge  of  the  wall.  In 
this  case  beam  theory  gives 

k - A D /h  (2) 

w 

In  the  second  case  (Fig.  2(b))  the  base  of  the  package  has  the  same 
flexural  stiffness  as  the  lid  and,  therefore,  with  the  package  under 
hydrostatic  pressure,  the  lower  edge  of  the  wall  will  be  subjected  to 
bending  moments  of  the  same  magnitude  as  the  upper  edge.  In  that  case 
It  can  be  shown  that 

k ««  2 D /h  (3) 

w 

The  actual  state  of  affairs  will  lie  somewhere  between  these  two,  and, 
therefore,  k will  be  given  the  following  representation: 

k “ a D /h  (A) 

w 
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where  a lies  between  4 and  2,  with  Its  actual  value  depending  on  the 
flexural  stiffness  of  the  base  of  the  package.  In  the  absence  of  better 
Information  a value  of  a equal  to  3 Is  suggested.  It  will  be  seen  later 
that  this  uncertainty  regarding  the  exact  value  of  a will  usually  be  of 
very  little  consequence. 

For  later  use  the  following  non-dlmenslonal  measure  of  k will  now 
be  Introduced: 


„ _ 4 ak  4.  w a. 

" = T 5 — 2 <“  r h> 

IT  C TT  C 


This  definition  can  be  rewritten  as  follows: 


(5) 


K H 


2 K 


(5a) 


where 


h 


c 

If  the  Poisson's  ratios  of  the  lid  and  wall  materials  are  the  same 
the  last  equation  reduces  to 


(6) 


w /Ws3  a 

“ r ¥ h 


(6a) 


c 

In  the  graphs  to  be  presented  it  will  be  advantageous  to  show  various 
quantities  as  functions  of  arctan  K,  rather  than  as  functions  of  K. 
Therefore,  Figure  3 Is  presented  to  permit  an  easy  conversion  from  K to 
arctan  K.  In  Figure  3 there  is  also  a graph  to  facilitate  the  conversion 
from  to  arctan  K. 
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BASIC  DATA  FROM  PLATE  THEORY 


Under  the  action  of  a uniform  gage  pressure  p (psi),  reactions 
will  develop  along  the  edges  and  at  the  corners  of  the  flat-pack  lid, 
as  depicted  in  Figure  h.  These  will  include  bending  moments  of  varying 
intensity  M (in. -lb.  per  in.)  due  to  the  restraint  against  rotation 
furnished  by  the  walla,  vertical  forces  of  varying  intensity  V(lb.  per  in.), 
and  concentrated  downward  forces  of  magnitude  R(lb.),  due  to  the  tendency 
of  the  corners  to  rise.  The  maximum  values  of  M and  V occur  at  the 
middle  of  the  long  side.  These,  together  with  R,  can  be  expressed  as 


M 

max 

" "l 

2 

• pa 

(7) 

V 

max 

• pa 

(8) 

R 

- n_ 

2 

• pa 

(9) 

where  n2  and  n^  are  functions  of  the  elastic  restraint  parameter  K 
and  the  aspect  ratio  b/a  of  the  lid.  The  values  of  n^^,  n2  and  n^ 
associated  with  any  given  configuration  can  be  obtained  from  Figures  5, 
6 and  7,  respectively.* 


*In  Figure  5 the  data  for  a clamped  plate  (K*“,  arctan  K » ir/2)  are  taken 
from  Table  35  of  reference  2.  All  other  data  in  this  figure  are  based  on 
the  analysis  in  appendix  A. 

In  Figure  6 the  data  for  a simply  supported  plate  (arctan  K ■ 0)  are  from 

Table  8 of  reference  2.  The  data  for  a clamped  plate  (arctan  K * it/2) 

are  based  on  the  analysis  in  appendix  B.  The  curves  for  elastically 
restrained  plates  (arctan  K •«  .4,  .8  and  1.2)  were  inserted  by  interpolation 
assuming  a linear  variation  of  n2  with  respect  to  arctan  K,  which  la  approx- 
imately the  variation  obtained  for  n^  and  n„.  In  view  of  the  small  change 
in  n2  in  going  from  simple  support  to  clamping  (around  6%  at  the  most)  and 
the  small  role  that  n„  will  play  in  the  subsequent  development,  this  linear 
interpolation  employed  in  Figure  6 is  considered  quite  accentable. 

In  Figure  7 the  data  for  a simply  supported  plate  (arctan  K = 0)  are  from 

Table  8 of  reference  2.  The  remaining  data  are  based  on  the  analysis  in 

appendix  A. 
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*The  data  in  Figure  8 for  the  cases  arctan  K = 0 and  tt/2  are  taken  froia 
Tables  8 and  35,  respectively,  of  reference  2.  The  remaining  data  are 
based  on  the  analysis  in  appendix  A. 

**Figure  9 applies  rigorously  only  to  the  case  K = 0,  v = 0.3,  but  ia 
probably  sufficiently  accurate  for  practical  purposes  ?or  other  values 
of  K and  v as  well.  The  solid  curves  of  this  figure  are  based  on  the 
table  on  page  181  of  reference  3.  The  dashed  curves  are  estimated 
interpolations  or  extrapolations  due  to  the  writer. 
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only  small  changes  In  arctan  K (see  Fig.  3)  and,  therefore,  only  small 
changes  In  n^,  n^,  n^  and  n^.  This  supports  the  assertion  made  earlier 
that  the  uncertainty  regarding  the  exact  value  of  a to  be  used  In 
equation  (A)  will  generally  not  have  serious  consequences. 


FORMULAS  FOR  TENSILE  STRESS  IN  THE  SEAL 


There  are  two  conceivable  locations  for  the  maximum  tensile  stress 
In  the  seal:  (a)  in  the  middle  of  the  long  side,  where  the  bending  moment 
transmitted  from  the  lid  to  the  wall  Is  a maximum;  and  (b)  at  the  corners, 
where  a reaction  force  R is  developed  to  prevent  the  corners  of  the  lid 
from  lifting  ur , We  shall  consider  these  two  locations  separately. 

Figure  10(a)  shows  the  assumedly  linear  stress  dlsliibutlon  across 
the  width  of  the  seal  at  the  middle  of  the  long  side  due  to  the  bending 
moment  M alone,  and  Figure  10(b)  shows  the  uniform  compressive  stress 

uuR  A 

at  the  same  location  due  to  Superposition  of  the  two  diagrams 


gives 


6 M V 

max  max 


for  the  maximum  tensile  stress  in  the  seal  at  the  middle  of  the  long 

side.  It  occurs  at  the  outside  edge  of  the  seal.  Substituting  for  M 

n 

and  V the  expressions  given  in  equations  (7)  and  (8),  one  obtains 
m&x 

s . “ p(-)^  n , 0 

edge  w odge 


where 


n j = 6n,  - n„  — 
edge  1 2a 


*Because  of  the  factor  w/a  in  equation  (13),  it  is  now  evident  that  ri2 
plays  only  a small  role  in  determining  the  magnitude  of  the  stress.  There- 
fore, the  simplifying  assumption  used  in  inserting  the  dashed  curves  in 
Figure  6 is  quite  acceptable. 
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*The  dimension  .47t  In  Figure  11  Is  based  on  some  rather  limited  theoretical 
data  from  p.  173  of  reference  2.  There  an  analysis  by  A.  Kromm  Is  cited, 
for  the  case  of  a simply  supported  square  plate  with  a/t  “ 20  and  with 
deflections  due  to  transverse  shear  taken  Into  account,  which  leads  to  a 
hlgh-lntenslty  distributed  reaction  In  place  of  the  concentrated  reaction 
R.  The  maximum  Intensity  of  this  distributed  reaction  In  that  particular 
example  can  be  correctly  obtained  by  assuming  the  concentrated  R of  clas- 
sical plate  theory  to  be  uniformly  distributed  of  a total  length  of  edge 
equal  to  .94  times  the  plate  thickness. 
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where 


n 5 n,/n, 

corner  3 6 


Equations  (12)  and  (19)  furnish  two  candidate  values  for  the 


maximum  tensile  stress  S in  the  seal.  Thus 

max 


S - p(-)^  n 
max  w max 


where 


n = the  larger  of  n . and  n 
max  edge  corner 


For  moat  designs  will  govern  in  equation  (22). 


APPLICATION  TO  DEfTGN  AND  SCREENING 


The  formulas  developed  above  can  be  of  use  both  to  the  designer, 

whose  objective  is  to  design  a package  that  will  remain  hermetic  under 

a specified  screening  pressure,  and  the  user,  whose  objective  is  to 

select  an  appropriate  screening  pressure  that  will  destroy  the  hermeticity 

of  packages  with  poor  quality  seals. 

Application  to  design;  Considering  first  the  designer,  let  it  be 

supposed  that  he  knows  the  sealing  material  to  be  used  and  has  a value 

of  the  allowable  tensile  stress  S for  that  material.  Then  his  criterion 

all 

for  a satisfactory  design,  from  the  point  of  view  of  retaining  hermeticity 
under  external  pressure,  should  be  that  as  given  by  equation  (21), 


is  less  than  S That  is, 

all  ’ 


p(-)^  n < S „ 
w max  — all 


where  p is  the  largest  anticipated  screening  pressure.  In  selecting 
the  designer  should,  of  course,  be  conservative.  If  is  taken  as  the 


median  tensile  strength  of  the  sealing  material • as  obtained  from  a number 
ot  tests  of  that  material  alone,  then  packages  designed  on  the  basis  of  the 
equality  sign  In  equation  (23)  will  have  a failure  rate  of  approximately 
50%  even  If  properly  sealed  (and  a higher  failure  rata  If  there  Is  a 
mixture  of  properly  and  Improperly  sealed  packages).  On  the  other  hand, 

If  the  designer  selects  for  the  1-percentlle  value  of  the  tensile 

strengths  obtained  In  a large  number  of  tents  of  the  sealing  material 
alone,  then  he  should  expect  only  a 1%  failure  rate  for  properly  sealed 
packages  designed  on  the  basis  of  equation  (23)  with  the  equality  sign. 

It  could  happen  that  the  designer  has  very  little  data  on  the 
distribution  of  tensile  strength  values  for  the  sealing  material,  or  even 
on  the  mean  strength,  but  he  doer  know  that  a certain  previously  designed 
package,  designated  as  I,  when  properly  sealed  with  the  same  material  had 
an  acceptable  failure  rate  F under  a screening  pressure  of  p^.  Then  In 
order  for  the  new  package,  designated  as  II,  to  have  a failure  rate  no 
greater  than  F when  properly  sealed  and  subjected  to  Its  screening  pressure 
p^_,  he  should  so  design  package  II  chat  Its  S Is  no  greater  than  that 

of  package  I.  Thus,  referring  to  equation  (21),  his  criterion  for  a 
satisfactory  design  of  package  II  should  be 


'II 


.a.  2 
(— ) n 
w max 


II 


i "I  ^ 


max  I 


(24) 


or 


/®\2  ^ ^I  /a\2 

(— ) n < (— ) n 

w max  — p w max 


(24a) 
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If  ^he  same  screening  pressure  is  to  be  used  for  package  II  as  for 
package  I,  this  design  criterion  becomes 


(241) 


Application  to  acreenltiR;  Turning  now  to  the  user  of  an  already 
designed  package,  let  us  first  suppose  that  he  has  a minimum  acceptable 
value,  S , for  the  tensile  strength  of  the  seal  and  wants  to  be 

sure  of  rejecting  all  packages  with  seal  strengths  less  than  that.  Then 
he  should  select  a screening  pressure  p such  that  the  S defined  by 

luaX 

equation  (21)  is  equal  to  or  greater  than  S . Thus,  the  required 

aCC0p  w 

screening  pressure  la  defined  by 


p > —acceet 

^ - /av2 


(25) 


—)  n 
w max 


It  Is  more  than  likely,  however,  c lat  the  user  of  the  package  will 
not  have  enough  Information  about  the  sealing  material  to  be  able  to 


specify  a value  for  S , but  he  might  know  that  in  the  past  a certain 

^ *'  accept 


pressure  p^  was  considered  suitable  for  the  screening  of  a certain  package, 
designated  as  I,  employing  the  same  sealing  material  as  the  present  package 
which  we  shall  designate  as  II.  Then  his  criterion  for  selecting  a suit- 
able screening  pressure  p^j.  for  the  present  package  should  b ' that  the 


of  package  II,  under  the  pressure  P-r^.,  be  at  least  as  big  as  the 


"max  ' ' * ‘II 

S of  package  I under  the  pressure  p . Referring  again  ..o  equation  (21), 


'max  ■ ‘ ■ ‘I 

we  see  that  the  criterion  for  selecting  p^^  is  then 


'II 


(26) 
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or 
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Possible  modification  In  above  procedures;  In  appendix  C some 
doubts  are  raised  as  to  the  validity  of  the  corner  stress.  Theoretical 
arguments  are  presented  which  Indicate  that  equation  (19)  might  over- 
estimate this  stress  and  that.  In  fact,  this  stress  might  not  exist  at 
all.  Thus,  there  Is  an  uncertainty  regarding  the  existence  of  the  corner 
stress,  and  the  question  naturally  arises  what  modifications.  If  any, 
are  needed  In  the  designing  and  screening  recommendations  just  presented 
because  of  this  uncertainty. 

In  the  writer's  opinion,  the  answer  to  this  question  should  be  as 

follows:  One  should  accept  or  discount  the  existence  of  the.  corner 

stress  on  the  basis  of  which  attitude  will  lead  to  the  more  conserva.tive 

result.  For  the  designer  this  means  that  if  he  is  designing  on  the 

basis  of  equation  (23)  he  should  use  that  equation  as  it  stands;  but 

if  designing  on  the  basis  of  equation  (24) , (24a)  or  (24b) , he  should 

make  sure  that  the  inequality  remains  satisfied  with  n replaced 

by  n , on  both  sides.  On  the  other  hand,  the  user  of  the  package, 
edge  

who  is  trying  to  arrive  at  a proper  screening  pressure  p,  would  make 
two  calculations,  one  on  the  assumption  that  the  corner  stress  exists, 
the  other  on  the  assumption  that  it  does  not  exist,  and  use  the  higher 
of  the  two  pressures  thus  obtained.  The  first  calculation  is  the 
one  already  described  above;  that  is,  it  is  a calculation  based  on 
equation  (25)  or  (26a) . For  the  second  calculation  equation  (25) 
or  (26a)  is  again  used,  but  with  replaced  by 


Fortunately,  for  most  packages  as  given  by  equation  (12), 

will  be  much  greater  than  given  by  equation  (19).  Thus,  in 

most  cases  the  question  as  to  whether  equation  (19)  over-estimates 

S or  whether  S actually  exists  will  be  academic,  and  the 

corner  corner  ^ ’ 

modified  designing  and  screening  procedures  will  be  in  no  way  different 
from  the  original  procedures. 

Lid  deflection:  It  goes  without  saying  that  the  designer  and  the 

user  both  must  concern  themselves  with  the  center  deflection  of  the  lid, 
in  order  to  Insure  that  during  any  screening  the  deflection  will  not  be 
so  large  as  to  produce  contact  between  the  lid  and  the  contents  of  the 
package.  The  pertinent  formula  for  estimating  center  deflection  is 
equation  (10a).  The  designer  should  so  design  the  package  that  under  the 
anticipated  screening  pressure  this  deflection  is  less  than  the  antici- 
pated internal  clearance  between  lid  and  contents.  By  the  same  token, 
the  user,  after  tentatively  selecting  a screening  pressure,  should  compute 
the  center  deflection  produced  by  that  pressure  and  make  sure  that  it  does 
not  exceed  the  internal  clearance.  If  it  does,  then  he  may  have  to  use  a 
less  severe  screening  pressure  (l.e. , accept  a higher  risk  of  non- 
hermeticlty)  in  order  not  to  jeopardize  the  integrity  of  the  contents  of 
the  package  in  the  screening. 

FLAT-PACKS  IN  A CENTRIFUGE 

As  part  of  the  total  screening  process,  packages  are  frequently  spun 
in  a centrifuge  in  such  a way  that  the  centrifugal  force  tends  to  push 
the  lid  into  the  cavity.  As  far  as  the  lid  alone  is  concerned,  this 
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centrifugal  force  Is  equivalent  to  a lateral  pressure  of  magnitude  Gtd, 

where  G Is  centripetal  acceleration  in  units  of  g (acceleration  of 

2 

gravity,  32.2  ft. /sec.  ) , t Is  the  lid  thickness,  and  d Is  the  npeclflc 
weight  (weight  per  unit  volume)  of  the  lid  material.  If  t Is  In  Inches 
and  d In  lbs.  per  cubic  Inch,  the  formula 

p - Gtd  (27) 

will  give  the  effective  pressure.  In  pal,  due  to  a centrifuge  acceleration 
of  G g's. 

Thus,  the  formulas  developed  In  the  preceding  sections  can  be  made 
to  apply  to  a package  In  a centrifuge  simply  by  replacing  the  symbol  p 
everywhere  by  Gtd.  In  this  way,  for  example,  the  following  formula  is 
obtained  from  equation  (10a)  for  the  central  deflection  of  the  lid  of  a 
package  in  a centrifuge: 


6 - 10.92  n,  ~ (f)^  • a • n^  (28) 

max  A E t 5 

c 

From  equations  (12) , (19)  and  (21)  the  following  formulas  are  obtained 
for  the  tensile  stresses  In  the  seal: 


S , - Gtd(-)^  n . 

edge  w edge 


S “ Gtd(— )^  n 

comer  w corner 


S ■“  Gtd(— )^  n 
max  w max 


(29) 

(30) 

(31) 


It  should  be  noted  that  the  Interaction  among  the  base,  walls,  and 
lid  of  a package  may  be  somewhat  different  for  a package  In  a centrifuge 
than  for  the  same  package  under  a uniform  hydrostatic  pressure.  There- 
fore, the  value  of  a to  be  used  in  equations  (4)  to  (la)  is  not  necessarily 
the  same  for  the  package  in  the  centrifuge  as  It  Is  for  the  same  package 
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under  hydrostatic  pressure.  As  an  aid  In  estimating  the  appropriate 
value  of  a let  us  consider  one  extreme  condition,  Illustrated  In  Figure 
14(a).  In  this  figure  the  lid  and  base  are  assumed  to  be  Identical  and 
the  package  Is  assumed  to  be  supported  In  the  centrifuge  only  at  Its 
edges.  Under  these  circumstances  the  lid  and  base  will  deform  Identically 
under  their  centrifugal  loadings  and  the  walls  will  develop  lines  of 
contraf lexure  at  mid-height . It  can  be  shown  that  the  appropriate  value 
for  a is  then  6.  A more  likely  support  condition  Is  illustrated  in 
Figure  14(b).  There  the  base  is  firmly  bonded  to  a flat  surface,  producing 
an  essentially  clamped  condition  for  the  bottom  edges  of  the  walls.  In 
this  case  a should  be  taken  as  4.  Again  it  should  be  noted  that  for  most 
packages  this  uncertainty  In  the  value  of  a is  of  no  real  consequence. 

Equations  (28)  to  (30)  can  be  used  to  determine  if  any  undesirable 
lid  deflection  or  seal  stresses  will  result  from  a given  centrifuge  accel- 
eration. Equation  (27)  can  be  used  in  the  following  way:  If  a certain 
hydrostatic  pressure  p is  known  to  be  destructive  of  lid  or  seal  for  a 
i Iven  package,  then  equation  (27),  rewritten  as 

C - ^ (32) 

will  give  the  acceleration  (in  g's)  that  will  produce  approximately  the 
same  destructive  effects.  Conversely,  If  a certain  hydrostatic  pressure 
is  known  to  be  non-destructive  of  lid  or  seal  then  equation  (32)  can  be 
used  to  estimate  a value  of  G that  will  also  be  non-destructive. 

To  gain  further  Insight  into  the  equivalence  of  hydrostatic  pressure 
and  constant  acceleration,  let  us  consider  a lid  of  .030  in.  thickness 
and  .302  lbs.  per  cu.  in.  density  and  ask  how  many  g's  of  centrifuge 


iPIlipippipgpffMVF^ 


acceleration  are  equivalent  to  a lateral  pressure  of  30  psi  on  the  lid. 
From  equation  (32)  the  answer  Is 

° " (.030) (.302)  “ 

NUMERICAL  EXAMPLES 


Example  1.  A package  has  the  following  properties:  a ■ b ••  1.000  in., 
t - .015  in.,  h - .105  in.,  w - .040  In.,  r - .010  in.,  E^  **  " 20,000,000 

3 

psl,  V "V  “O.l.d*  .302  lb. /in.  . What  are  the  maximum  tensile  stress 
c w 

in  the  seal  and  the  central  deflection  of  the  lid  due  to  an  external  pressure 
of  30  psl? 

With  Cl  assumed  to  be  3,  equation  (6a)  gives 


K 


1 


’.040\^  1.000 
.015  .105 


- 541 


(33) 


Figure  3 then  gives  arctan  K •“  1.57  (that  Is,  the  lid  is  essentially 
clamped  along  its  edges).  Enten.-t;  Figures  5,  6,  7 and  8 with  b/a  » 1 
and  arctan  K <»  1.57,  one  obtains 

» .051,  n^  “ .443,  n^  “ “ .00125  (34) 


For  p >■  30  psl. 


4 

_£2_ 

4 

E t 
c 


30 

20,000,000 


29.6 


Therefore  (Fig.  9) 

n_  - .825 

i 

Equation  (10a)  no\?  gives  the  following  lid  deflection: 


(35) 


(36) 


6 

max 


10. 92(. 00125) 


30 

20,000,000 


(1.000  ln.)(.825  ) 


“ . 0050  in. 


(37) 
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From  equations  (13),  (17a),  (20)  and  (22)  there  follows 


‘^edge  ’ 1?^ 


.288 


1 + 


- •-  (M 


1.34 


(38) 


corner 


max 


.288 


Thus,  the  maximum  tensile  stress  in  the  seal  occurs  at  the  middle  of  the 
long  edge.  Its  value,  from  equation  (21),  is 

2 


(39) 


Example  2,  Suppose  that  30  psl  is  considered  a satisfactory  screening 
pressure  for  the  package  of  Example  1 (to  be  referred  to  as  package  I), 
what  screening  pressure  would  be  appropriate  for  a second  package 
(package  II)  identical  in  all  respects  to  package  I except  for  the  length 
b,  which  is  2.000  in.?  What  will  the  lid  deflection  be  under  this  screen- 
ing pressure? 

For  package  II  we  have  b/a  •=  2 and  arctan  K still  equal  to  1.57. 
Figures  5,  6,  7 and  8 then  give  the  following  Information  for  this  package: 


nj^  ■ .083,  n2  ■ .51,  n^  ^ 0, 


.00255 


(40) 


A1l)o,  from  equations  (13),  (17a),  (20)  and  (22), 


"edge  - riso  - 


"6  ■ 


(41) 


n 


corner 

= .477 


max 
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Equation  (26a)  now  gives 


PlI  " 


. 2 

fl.000\  . 

rr^j 


■ 18  psi 


as  the  appropriate  screening  pressure  for  package  II.  For  this  pressure 


18  /l.000\^  , 

J.4  " 20,000,000  V *015/  “ 


The  curve  for  b/a  - 2 in  Figure  9 then  gives 

n^  - .87  (AA) 

From  equation  (10a)  the  center  deflection  of  the  lid  of  package  II  under 
a screening  pressure  of  18  psi  is  now  found  to  be 


'max  ■ loTosboO 

= .0065  in. 


(1.000  ln.)(.87) 


Example  3.  For  the  package  of  Example  1 a maximum  tensile  stress 

of  5,A00  psi  in  the  seal  material  was  computed  for  a screening  pressure 

of  30  psi.  Suppose  that  the  maximum  permissible  stress  in  the  seal  is 

A, 500  psi,  what  new  lid  thickness  (if  any)  would  reduce  S to  this 

value,  all  other  dimensions  remaining  unchanged? 

Assume,  tentatively,  that  will  still  occur  along  the  edge, 

rather  than  at  the  corner.  Than  the  criterion  for  a satisfactory  design, 

from  equation  (23)  with  n set  equal  to  n , , is 

^ max  ^ edge 


edge  — 2 

,a. 


A500 
/1.000\2 
*(  .oAoy 


In  view  of  equation  (13)  this  criterion  may  be  written  as 


6n,  - n-  — < .24  (45) 

12a  — 

or,  substituting  for  w/a  Its  value  of  .040  and  solving  for  n^, 

nj^  £ .040  4-  .0067  n^  (46) 

The  value  of  Is  fairly  Insensitive  to  K (see  Fig.  6)  and,  therefore, 
very  little  error  is  Involved  If  we  continue  to  use  for  the  value 
found  for  It  In  Example  1,  namely  .443.  Thus,  equation  (46)  becomes 

n^  £ .043  (47) 

Referring  to  the  curve  for  b/a  =■  1 In  Figure  5,  we  see  that  this  requires 
arctan  K < 1.45  (48) 


or,  from  Figure  3, 


K < 20 


(49) 


Substituting  for  K Its  expression  from  equation  (6a),  and  using  the 
numerical  data  specified  in  Example  1,  we  convert  this  to 


whence 


or 


3(f)  20 


f 1 .888 


(50) 


(51) 


^ ^ w .040  . 

^ - .888  ’■  .888  “ 


Thus,  a lid  thickness  comparable  to  the  wall  thickness  would  be 
needed  to  reduce  the  maximum  tensile  stress  in  the  seal  at  the  middle 
of  the  side  to  4500  psi. 

The  corner  stress  in  the  seal  must  now  be  checked  to  see  that  It 
does  not  exceed  the  maximum  permissible  stress  of  4500  psi.  From 
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equations  (19),  (20)  and  (17a),  we  have 


S “ p (— ) ^ 

corner  ^ w 


(52) 


1 + .94^-  .21(1) 

For  arctan  K “ 1.45  and  b/a  - 1.0,  Figure  7 gives 


n^  “ .0056 


(53) 


Substituting  this  into  equation  (52)  together  with  the  numerical  values  of 
the  other  parameters,  we  have 


S 


corner 


30 


/l.000\^ 1 

V.OAoj 


0056 


097) 

0401 


l\.  2l(^\ 
oj  • \.04oj 

which  is  Indeed  less  than  the  4500  psl  allowable. 


2“  - 32  psl 


(54) 


CONCLUDING  REMARKS 


Formulas  based  on  classical  plate  theory  have  been  presented  for  the 
reactions  and  center  deflection  of  the  lid  of  a rectangular  flat-pack 
under  external  pressure,  and  suggestions  have  been  made  as  to  how  these 
formulas  can  be  of  use  in  the  design  and  screening  of  flat-packs. 

There  are  several  simplifying  assumptions  underlying  the  present 
work.  Most  of  them  have  already  been  discussed  and  justified.  One 
which  has  not  yet  been  discussed,  and  which  warrants  close  scrutiny,  is 
the  second  hypothesis  stated  in  the  Introduction;  namely,  that  the  screen- 
ing effectiveness  of  external  pressure  in  weeding  out  non-hermetlc  or 
potentially  non-hermetic  seals  is  due  to  the  tensile  stresses  which  this 
pressure  produces  in  the  seal  material.  This  hypothesis  is  admittedly 
debatable.  One  could  argue  that  the  screening  action  is  due  to  the 
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seepage  of  fluid  into  the  package  through  pre-existing  voids  or  narrow 

passages  in  seals  of  poor  quality.  Undoubtedly,  some  of  the  screening 

action  is  due  to  such  seepage,  but  the  writer  feels  that  the  tensile 

stresses  produced  by  the  pressure  are  the  predominant  screening  mechanism. 

There  is  some  evidence  to  support  this  feeling:  Firstly,  those  who  are 

knowledgeable  in  the  screening  of  flat-packs  report  that  the  screening 

power  (failure  rate)  of  external  pressure  goes  up  drastically  as  the 

package  size  (dimensions  a and  b)  increases.  This  is  consistent  with 

equation  (21) , which  shows  that  the  maximum  tensile  stress  in  the  seal 

is  approximately  proportional  to  the  square  of  the  dimension  a.  Secondly, 

p.  5 of  reference  4 describes  the  "pressure  withstanding  capability"  of 

properly  sealed  Bendix  packages  as  follows: 

Holder  Size  Pressure  Differential 

3/8  square  5 atmospheres 

3/8  to  5/8  square  3 atmospheres 

5/8  to  3/4  square  2 atmospheres 

3/4  square  1.25  atmospheres 

We  note  that  as  the  size  is  doubled  (going  from  3/8  square  to  3/4  square), 

the  allowable  pressure  differential  is  reduced  by  a factor  of  4 (from  5 

atmospheres  to  1.25  atmospheres).  Equation  (26a)  predicts  approximately 

the  same  relationship  between  package  size  and  screening  pressure. 

A second  assumption  which  warrants  further  attention  is  the  one 
illustrated  in  Figure  10(a);  namely,  that  the  bending  moment  at  the  edge 
of  the  lid  produces  a linear  variation  of  stress  across  the  width  of  the 
seal,  as  in  elementary  beam  theory.  This  kind  of  assumption  is  frequently 
made;  it  follows  from  the  assumption  of  plane  horizontal  sections  through 
the  seal  remaining  plane  and  the  assumption  of  Hooke's  law.  As  the 
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stresses  in  the  seal  material  becomes  high  we  can  expect  that  plastic 
behavior  and  possibly  creep  in  the  seal  material  will  destroy  the 
linearity  in  the  stress  distribution  and  reduce  the  maximum  stress. 
Furthermore,  if  the  lid  is  thin  compared  to  the  thickness  of  the  wall 
(t<<w),  we  should  not  expect  the  lid  to  be  able  to  enforce  the  plane 
sections  assumption  in  the  seal.  For  t<<w  the  action  of  that  portion 
of  the  lid  in  contact  with  the  seal  would  probably  be  like  that  of  a 
beam  attached  to  an  elastic  foundation,  ana  the  pressure  distribution 
between  such  a beam  and  its  foundation  is  known  to  be  non-linear.  Thus, 
the  linear-stress-distribution  assumption  of  Figure  10(a)  has  obvious 
shortcomings.  Nevertheless,  it  has  been  employed  in  the  hope  that  it 
may  be  sufficiently  good  for  practical  purposes  and  because  any  more 
sophisticated  assumption  would  involve  an  inordinate  amount  of  additional 
analysis . 

A third  assumption  implicit  in  the  present  work  is  that  the  pressure 
on  the  side  walls  makes  a negligible  contribution  to  the  stresses  in  the 
seal.  For  shallow  packages  the  validity  of  this  assumption  is  intuitively 
evident.  A semi-quantitative  justification  of  it  can  be  made  by  analyz- 
ing the  configuration  of  Figure  1(b)  as  a rigid-jointed  plane  frame  under 
pressure  on  the  vertical  members  alone  and  pressure  on  the  horizontal 
members  alone.  Such  an  analysis  shows  that  the  ratio  of  the  corner  bending 

moment  due  to  side-wall  press’'‘..e  alone  to  the  corner  bending  moment  due  to 

3 3 

pressure  on  the  top  and  bottom  covers  alone  is  of  the  order  of  D^h 
3 3 3 3 

or  (E^/E^)  (t  /w  ) (h  /a  ) . For  any  flat-pack  of  practical  design  this 
ratio  will  be  much  smaller  than  unity,  indicating  that  the  side-wall 
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pressure  makes  a negligible  contribution  to  the  bending  moment  transmitted 
across  the  seal  and,  therefore,  to  the  stress  in  the  seal. 

The  final  test  of  the  validity  of  any  theoretical  work  and  the 
hypotheses  on  which  It  is  based  must,  of  course,  be  based  on  experiments. 
In  the  present  case  one  can  envision  a modest  experimental  program  that 
could  provide  evidence  as  to  the  validity  of  the  assumptions  that  have 
been  made.  This  program  would  require  flat-packs  of  different  sizes 
and  shapes,  all  carefully  sealed  with  the  same  material.  The  value  of 
n (eq.  (22))  would  be  computed  for  each  package.  Each  package  would 
then  be  subjected  to  external  pressure  which  is  Increased  in  small  steps, 
and  the  pressure  p^^  which  first  causes  a loss  of  hermetlcity  would  be 
noted.  The  product  p^^  (a/w)  would  then  be  computed  for  each 

package.  According  to  equation  (21),  this  product  represents  the  maximum 
tensile  stress  in  the  seal  at  the  time  of  loss  of  hermetlcity.  If  the 
seals  are  of  reasonably  uniform  quality,  so  that  they  all  have  essentially 
the  same  tensile  strength  and  if  the  hypothesis  of  this  report  (to 

the  effect  that  fai].ure  occurs  when  S ■ S - ) is  correct,  then  this 
product  should  turn  out  to  be  approximately  the  same  for  all  packages. 

If  the  test  program  Just  described  tends  to  confirm  the  present 

hypotheses,  then  a second  program  would  be  desirable,  identical  to  the 

first,  except  that  the  seal  material  would  be  changed  to  one  of  much 

higher  or  much  lower  strength.  The  purpose  of  this  program  would  be  to 

see  if  there  is  a corresponding  change  in  the  experimental  values  of 
2 

p (a/w)  n , as  there  should  be  if  the  failure  hypothesis  S “ S - 
*^cr  max’  max  ult 

is  correct. 
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APPENDIX  A 


SMALL-DEFLECTION  ANALYSIS  OF  A UNIFORMLY  LOADED  ELASTIC 
RECTANGULAR  PLATE  WITH  EDGES  ELASTICALLY  RESTRAINED  AGAINST  ROTATION 


Here  the  energy  method  will  be  employed,  in  conjunction  with 
double  Fourier  series  and  Lagranglan  multipliers,  to  obtain  the 
deflections  and  other  quantities  of  Interest  in  a uniformly  loaded 
elastic  rectangular  plate  with  edges  elastically  restrained  against 
rotation. 

Coordinate  system  and  notation.-  We  shall  employ  the  coordinate 
system  and  notation  shown  in  Figure  4 and  let  w(x,y)  denote  the 
deflection  of  the  plate,  positive  if  in  the  same  direction  as  the 
pressure,  p.  The  following  additional,  symbols  will  be  used: 

k = elastic  restraint  stiffness  (moment  per  unit  length  per 
radian -of  rotation). 

E - Young's  modulus  of  plate  material. 

V = Poisson's  ratio. 

D = Et^/[12(l-v^)] . 

K = 4ak/(7t^D) . 

B s a/b. 

Series  for  w and  its  second  derivatives.-  The  double  Fourier  series 
M N 

w(x,y)  = y y a sin(mTTx/a)  sin(mry/b)  (Al) 

tD*  X 1 3 ^ t * 1 } 3 ^ • • 

is  an  appropriate  representation  of  the  deflections  of  a rectangular 
plate  with  zero  deflection  at  the  edges.  Only  odd  values  of  the 


-24- 


summation  Indices  m and  n are  taken  In  equation  (Al)  because  under 
uniform  pressure  the  deflections  will  be  symmetric  with  respect  to 
the  planes  x-a/2  and  y-b/2.  Upper  limits  M and  N are  shown  for 
the  summations  with  the  understanding  that  In  the  computations  M and 
N will  be  assigned  sufficiently  high  values  to  establish  practical 
convergence  for  all  quantities  of  Interest. 

The  boundary  condition  of  zero  deflection  at  the  edges  permits 
2 2 2 2 2 

one  to  obtain  3 w/3x  , 3 w/3y  , and  3 w/3x3y  by  termwise  differentiations 
of  aquation  (Al) . Consequently, 


3x3y 


M 

I 

m-1,3, . . 
M 

I 

m*l,3, . . 

M 

I 

m*l , 3 , . . 


2 

y a (mir/a)  sln(mTix/a)  sln(mry/b) 

n"l,3, . . 


• 2 

y a^(mt/b)  sln(mTix/a)  8ln(mry/b)  (A2) 

n*l , 3 , . • 


N 

y a (mTt/a)  (mr/b)  cosCmirx/a)  cos(mry/b) 

‘'o  mn 


Series  for  edge  rotations.-  The  rotations  Oj^(x)  of  the  edge  y*0 


and  02(y)  of  the  edge  x«0  will  be  represented  In  the  form  of  single  series 


0^(x) 


M 

) b sin(nnrx/a) 

m-1,3,..  “ 


02^y) 


N 

y c ain(miy/b)  (A3) 

n-1,3, . . 


Constraining  relations.-  The  deflections  w(x,y)  and  the  rotations 
0^(x)  and  02 (y)  must  satisfy  the  following  compatibility  conditions: 

(3w/3y)y.O  “ (9w/3'0^^q  “ 02(y) 
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These  give  rise  to  the  following  constraining  relations  among  the 

Fourier  coefficients  a , b , c : 

mn  m’  n 

^ - 0 (m-1,3 M) 

mn  b m 

(A4) 


I 

n-1,3 


M 


I 

m“l , 3 , . . 


miT 

a — 
mn  a 


c 

n 


0 (n-l,3,...,N) 


Strain  energy  of  plate.-  The  strain  energy  Integral  for  a 
rectangular  elastic  plate  undergoing  small  deflections  Is  given  on 
page  342  of  reference  2.  There  It  is  also  shown  that  substitution  of 
the  series  (Al)  into  the  strain  energy  Integral  leads  to  the  following 
formula  for  the  strain  energy  of  the  plate: 


TT^abP  y y 2 , m^  n^.2 

® m-1,3,..  n-1,3,..  b^ 


(A5) 


Strain  energy  of  elastic  restraint.-  The  strain  energy  V2  of  the 
rotational  elastic  restraint  along  the  edges  of  the  plate  can  be  written 


as 


V2  - 21  y k 


[0^(x) dx  + ^ k 


Substitution  of  the  series  (A3)  for  0j^(x)  and  62 ^y^  leads  to 


2 2 


, M N 

- ka  I b^  + y kb  I 


(A6) 


1r%  m *d  10 

, 3 , . . n-1 , 3 , . . 

Potential  energy  of  the  applied  load.-  Taking  the  xy-plane  taken  as 
datum,  one  may  express  the  potential  energy  of  the  applied  load  as 

- -p  w dx  dy 
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I. Ii|i|ivmnp  I 


or,  substituting  for  w its  series  expression. 


4pab 


M N a 

V V tnn 

1-3  1 O 


(A7) 


Minimization  of  the  total  potential  energy.-  The  total  potential 
energy  (TPE)  of  the  system  consisting  of  the  plate,  its  loading,  and  the 
elastic  constraints  is  defined  as 


TPE  =■  + V2  + 


(A8) 


According  to  the  principle  of  minimum  total  potential  energy  for  elastic 

structures,  the  true  configuration  of  the  system,  among  all  possible 

configurations  satisfying  continuity  and  any  rigid  constraints,  is  that 

which  minimizes  the  TPE.  Thus,  the  correct  values  of  the  a , b and  c 

mn  m n 

are  those  which  minimize  expression  (A8)  while  at  the  same  time  satisfying 
the  constraining  relations  (A4).  This  "constrained  minimization"  problem 
can  be  solved  by  setting  up  the  new  function 


TPE*  = V,  +V-+V„-  y ct(  I a ^-b) 
1 2 3 - m , mn  b m 

ni*X  f 3 ^ 1 3 9 • • 


N 

I 

n=l,3, . 


M 

I 


miT 


mn  a 


c_) 


(A9) 


m-1,3, . 

where  the  and  3^  are  Lagranglan  multipliers,  and  making  it  stationary 

with  respect  to  the  a 's,  the  b 's,  the  c 's,  the  u 's,  and  the  3 's 

mn  m n m n 


via  the  following  conditions: 

3(TPE*)/aa^j  - 0 (i-1,3 M;  J-1,3,...,N) 

a(TPE*)/3b^  - 0 (i-l,3,...,M) 

3(TPE*)/3Cj  - 0 (j-l,3,...,N) 


(AlO) 

(All) 

(A12) 


I 
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3(TPE  )/3a^  - 0 
3(TPE*)/33j  “ 0 


(1-1, 3,..., M) 

(j-1,3 N) 


(A13) 

(AlA) 


(Equations  (A13)  and  (A14)  merely  replicate  the  constraining  relations 
(A4).)  With  the  differentiations  carried  out  and  some  minor  algebraic 


manipulations  executed,  (AlO)  to  (A14)  become: 

(1— 1,3,  • • • ,Mj  j— 1,3,  •••!)) 


IJ 

* 

*1  ■ 
it 

■ 

'j 


2 2 
1 + + 1“^^ 


ij(l^+J^6^)^ 


(A15) 


- Aj^/K 

(1-1, 3,..., M) 

(A16) 

- V“ 

(J-1,3,...,N) 

(A17) 

19  I 

J— 1, 3 , • 

JAij  - 0 

• 

(1-1, 3,..., M) 

(A18) 

M 

j 1 

1-1,3,.. 

- 0 

(j— 1,3, • » • ,N) 

(A19) 

V it  if 

where  A^j,  A^,  , b^,  Cj,  and  K are  dimensionless  measures  of  a^^ , a^, 

b^,  Cj,  and  k,  defined  as  follows: 


IJ 


m ti^Da^^j/lba^p 

(A20) 

m lTi^a^/4ab^p 

(A21) 

= jTi^8j/4a^bp 

(A22) 

y ln^Db^/16ab^p 

(A23) 

- jTT^DCj/16a\ 

(A24) 

2 

! 4ka/n  D 

(A25) 
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Equations  (A15)  through  (A19)  are  a system  of  simultaneous  linear 

A it 

algebraic  equations  defining  the  , A^,  Bj , b^,  and  Cj . 

Simplification  of  simultaneous  equations.-  Equations  (A15),  (A16) 

it  it 

and  (A17)  can  be  used  to  eliminate  the  A^^ , b^  and  Cj  In  (A18)  and  (A19) 
The  latter  then  becomes: 


2 2 
1 + j‘‘A^  + 1^‘b 


^1  + 3K  \ 2222 


0 (1-1, 3,... M)  (A26) 


M 1 + i\  + ±h 

B + K I T2tl ■ ° 

J 1-1,3,..  (l^+j'^3^)^ 


(J-1,3,...,N)  (A27) 


A further  reduction  In  the  number  of  simultaneous  equations  Is  effected 


as  follows:  Solve  (A27)  for  each  Bj  in  terms  of  all  the  A^  to  get 

f2<J)  ^2  M A 

“j  ■ fi(J)  ■ t^(J)  (,2+jV)'  


where 


f (j)  m k"^  + i 

1-1,3,.. 


(A29) 


f,(j)  = - i (i^+j^3^*)"^ 

1-1,3,.. 

Now  use  (A28)  to  eliminate  the  B^'s  in  (A26)  and  thereby  arrive  at  the 
following  system  of  simultaneous  equations  in  the  A^'s  alone: 


I A f^(i,q)  - f^(i) 


(1-1,3 M) 


• (A30) 


q-1,3,. . 


where. 


N 

I 


2 2 22-2 


- I (j)]"^  (A31) 

J-1,3,.. 


f.d)  = - I {l+i^f,(J)[fi(J)]“-^}  (A32) 

J-1,3,.. 

with  6^^  being  Kronecker's  delta. 

On  the  basis  of  these  tr  avis  format  ions,  we  now  have  the  following 
calculation  procedure:  Solve  equations  (A30)  (whose  number  depends  only 

on  M,  not  on  N)  simultaneously  for  the  substitute  their  values 

into  (A28)  to  obtain  the  B.'s.  (In  the  special  case  k*0  the  A 's  and 

J q 

the  B 's  are  identically  zero  (see  eqs.  (A26)  and  (A27)).  With  the  a 's 

J q 

and  Bj's  known,  equations  (A15)  will  yield  the  values  of  the  A^j's. 

Once  the  above  quantities  have  been  determined,  all  other  items  of 
interest  can  be  computed,  as  will  be  seen  in  the  following  sections. 

Center  deflection.-  Equation  (Al) , evaluated  at  x*a/2,  y=b/2,  gives 

m+-n-^ 

(A33) 


,-li 


N 


«(f.  |)  - ! I (-1) 

i * Ot>  1 o in  nm 

TT  E ra=l,3,..  n=l,3,.. 


Corner  reaction.-  From  page  85  of  reference  2 we  have  the  following 
equation  for  the  corner  reaction  in  terms  of  the  corner  rate  of  twist: 


R = 2D(l-v) (3  w/3x3y) 


x=0 

y»0 
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Elimination  of  the  twist  through  equation  (A2)  then  gives 


R » pa^«  32(1-v)Ett  ^ \ I A mn 


M 


N 


(A34) 


Bending  moments  along  the  sides.-  The  bending  moments  normal  to  the 
edges  y-0  and  x"0  are  related  to  the  rotation  functions  6j^(x)  and 
as  follows: 


(M  ) - -kO  (x) 

^ y-0 


(A35) 


(M  ) - -ke„(y) 

* x-Q  ^ 


Substituting  for  0j^(x)  and  their  series  expansions  (A3)  and  making 

use  of  the  relations  (A16)  and  (A17),  one  obtains 


(M  ) . 

^ yV“0  i . , 

IT  m-1,3,.. 


M 

I 


A 8in(mnx/a) 
m 


(A36) 


. 2 N 

(M  ) „ - ■ - 'l 

X x-0  3 

IT  n-1,3,. 


sin(niTy/b) 


Curvatures  along  the  sides.-  The  curvatures  normal  to  the  edges 
can  now  be  determined  from  the  moment-curvature  relations  and  equations 
(A36)  as  follows: 


( - i(M  ) - - ^ 

' 2'  d'  v'v=0  3 

9y  y=0  ^ ^ IT  u m-1,3,.. 


M 

I 


A sin(mTTx/a) 



m 


(A37) 


-4<Vx-0 

3x  x-0 


- A£§_ 


N 

I. 


II  D n-1,3,, 


B^  8in(niiy/b) 
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Some  third  derivatives  of  w.-  Termwise  differentiation  of  the 


Fourier  series  (Al)  is  permissible  in  order  to  arrive  at  the  following 
series  expansions: 


Dx3y  m*l,3,..  n-1,3,.. 


/nil,,  2, mil.  miTx  , n^W 

a (—r)  ( — ) cos  sin  —r^ 

mn  b a a b 


(A38) 


T " I ^ “<>»  ^ 

f\  In  ^ n I&n  3D  3 D 


3 3 

Termwise  differentiation  is,  however,  not  valid  for  obtaining  9 w/9x 
3 3 

and  3 w/3y  . We  therefore  postulate  the  new  series 


3 M N 

3 w V r 1.  niiix  , mi 

— r » / A b cos  sin  — ^ 

.„3  mn  a b 

3x  m*‘l,3, . . n“l,3, . . 


.3  M N 

I I 

3y  m-1,3, . . n»l, 3, . 


V , mnx  mi 

A c sin  cos 

*-*  mn  a 


(A39) 


from  which  it  follows  that 


4 3 w mux  , miy  . , 

3 ” — — ;r  cos  sin  — r-^  dx  dy 

mn  ab  ,,  3 a b 

0 0 


(A40) 


"mn  ab 


3 w miy  . mux  . , 

— r cos  — sin  dy  dx 

,,  3 b a ^ 


/a'.  <■/■■!- i’xi. 


Partial  Integration  with  respect  to  x In  the  first  of  these  equations 
gives 


b . _ A [ 2 

mn  a b 


(■^^)  sin  dy] 
3x  x-0 


+ — I ^ 

a ab 


3 w . imrx  . niry  . , , 

— _ sin  sin  — ^ dx  dy] 

r\  D 

3x 


(A41) 


The  first  bracketed  term  will  be  recognized  as  the  formula  for  the 

2 2 

Fourier  coefficient  In  the  sine  series  for  (3  w/3x  ) and  therefore, 

x"0 

referring  to  equations  (A37)  , it  may  be  replaced  by 


4pa^  ^ 

■ IT^D  * " 

The  second  bracketed  term  Is  the  formula  for  the  Fourier  coefficient 

2 2 

In  the  slne-slne  series  expansion  for  3 w/3x  , and,  referring  to 

2 

equations  (A2) , It  may  be  replaced  by  -a^(mTr/a)  . As  a result,  equation 
(A41)  becomes 


, 16pa  n /mif\3 

“ mn  a 


mn  3„  n 

ir  D 


(A42) 


Analogous  operations  on  the  second  of  equations  (A40)  give 


16pb  m 


mn 


tt^D 


,mi.3 

m ^mn  b 


(A43) 


Thus,  the  b and  c in  equations  (A39)  are  now  expressed  in  terms 
mn  mn 

of  quantities  which  are  known  once  the  A^j , A^  and  have  been  computed. 
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Computational  details.-  The  above  analysis  was  the  basis  for 
computing  the  and  n^  data  in  Figures  5 and  7 for  values  of  K other 
than  0 and  Calculations  were  made  for  M « N ■ 11,  13,  15  and  the 
results  were  graphically  extrapolated  to  M = N » «>.  This  technique 
could  not  be  used  successfully  for  computing  n2  because  of  poor 
convergence,  or,  perhaps,  non-convergence,  of  the  series  (A45) . 
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APPENDIX  B 


SMALL-DEFLECTION  ANALYSIS  OF  A UNIFORMLY  LOADED  ELASTIC 
RECTANGULAR  PLATE  WITH  CLAMPED  EDGES 

Here  the  energy  method  is  again  employed  in  conjunction  with 
Lagrangian  multipliers  but  with  a different  double  Fourier  series,  to 
analyze  a uniformly  loaded  rectangular  plate  with  clamped  edges.  The 
main  purpose  is  to  obtain  data  for  the  K = » curve  of  Figure  6 (which 
I data  could  not  be  obtained  from  the  analysis  in  appendix  A because  of 

poor  convergence).  At  the  same  time,  expressions  are  given  for  the 
deflections  and  the  edge  bending  moments.  It  will  be  seen  that  the 
data  needed  for  the  K = “>  curve  of  Figure  6 are  obtained  directly 
from  the  values  of  the  Lagrangian  multipliers. 

Coordinate  system  and  notation*-  The  same  coordinate  system  is 
used  as  in  appendix  A,  and  the  symbols  k,  E,  v,  D,  K and  3 will  have 
the  same  definitions  as  in  appendix  A. 

Series  for  w.-  The  following  series  gives  zero  normal  slope  at 
the  edges  and  the  deflection  symmetry  appropriate  to  a uniform  loading; 

M N 

w(x,  y)  = J 'I  a cos(mTTx/a)cos(mry/b)  (Bl) 

m=0,2,. . n=0,2,. . 

It  i.s  understood  that  calculations  will  be  made  with  successively  higher 
values  of  M and  N until  convergence  has  been  achieved  for  all  physical 
quantities  of  interest. 

Constraining  relations.-  The  conditions  of  zero  deflection  along 
the  edges  requires  that  the  following  constraining  relations  be  imposed 
on  the  coefficients: 


0 (n  * 0,  2,  4,  , 


N) 


(B2) 


M 

I 

m»0 ,2, . 


mn 


N 

y a "0  (m»*0,2,4,...,M)  (B3) 

no 

n*0 , 2 , . . 


These  equations  are  redundant,  inasmuch  as  summing  (B2)  over  n leads 
to  the  same  equation  as  summing  (B3)  over  m.  In  order  to  remove  this 
redundancy,  we  shall  replace  the  first  of  equations  (B2)  and  the  first 
of  equations  (B3)  by  the  single  equation  obtained  by  summing  these  two. 
Thus,  the  following  set  of  independent  constraining  relations  is 
obtained  to  replace  (B2)  and  (B3) : 


2a 


N 


00  ««7  A ^On  ^ , “mO 

n*2,4, . . m“2,4, . . 


M 

I 


- 0 


M 

I “ 0 (n  = 2,  4,  6,  ...,  N) 

^ n lun 

m*0. 2, , . 


N 

I ^Tnr.  “ 0 (m  - 2,  4,  6,  ...,  M) 
n==0,2, . . 


(B4) 


(B5) 


(B6) 


Strain  energy  of  plate.-  Substitution  of  the  series  (Bl)  into 
the  strain  energy  expression  for  a rectangular  elastic  plate  undergoing 
small  deflections  (p.  342  of  ref.  2)  gives  the  following  formula  for 
the  strain  energy  of  the  plate: 


I 


8a  m=0,2,..  n=0,2,.. 


I (»"  + + «.o  + 


where  6 . and  are  Kronecker  deltas. 
mO  On 
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Potential  energy  of  the  applied  load*-  With  the  xy-plane  as 


datum,  the  potential  energy  of  the  applied  load  Is 
■“  -p  //  wdxdy 

or,  in  view  of  equation  (Bl), 

Minimization  of  the  total  potential  energy,-  The  total  potential 
energy  (TPE)  is  defined  as 

TPE  “ (B9) 

By  the  same  argument  as  in  appendix  A,  the  correct  values  of  the 
in  the  series  (Bl)  are  those  which  minimize  expression  (B9)  while  at 
the  same  time  satisfying  the  constraining  relations  (B4)  to  (B6) . 

Again  we  have  a "constrained  minimization"  problem  which  can  be  solves’ 
by  setting  up  the  new  function 


TPE.  - V,  + V3  - y(2a  + f 

n=2,4,..  m=“2,4,.. 


N M M N 

~1  a)“I  3) 

n=2,4,.."m=0,2,.  m-2,4,..“  n=0,2,..™" 


(BIO) 


where  Y,  are  Lagrangian  multipliers,  and  making  it  stationary 


with  respect  to  the  a 's,  Y,  the  a 's,  and  the  3 's  via  the  conditions 

mn  n m 


3(TPE*)/3aQQ  = 0 
3(TPE*)/3a^^  - 0 
3(TPE*)/3aQj  » 0 
3(TPE*)/3a^j  « 0 


(i  = 2,  4,  6,  ...,  M) 

(j  “ 2,  4,  6,  . . . , N) 

(i  •=  2,  4,  ...,  M;  j « 2,  4, 


(BID 
(B12) 
(B13) 
N)  (BL4) 
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3(TPE*)/9y  ■ 0 


(B15) 


3(TPE*)/8S, 


- 0 

(.1  - 2,  4,  6, 

....  N) 

(B16) 

- 0 

(1  - 2,  4,  6, 

....  M) 

(B17) 

The  last  three  of  these  conditions  merely  reproduce  equations 

(B4) , (BS)  and  (B6) . With  the  differentiations  carried  out,  the  first 

four  conditions  read  as  follows: 


-pab  - 2y  0 


(B18) 


TT^D  . 1^  . - « 

8ab  2 ^^10  ■ Y ~ 

P 


0 (1  " 2,  4 M) 


(B19) 


A . ,4  2 
Sab  J ^ 


4a 


Oj 


Y - otj  - 0 (j  - 2,  4 N)  (B20) 


ir^D  1^  2 2 “ 2,  4,  . . . , M. 

8^  ("y  + j 8)  '2a^j  -Oj-B^-0  j-2,4,  ...,N 


(B21) 


Equations  (B18)  through  (B21) , together  with  (B4),  (B5)  and  (B6),  are  a 
system  of  simultaneous  linear  algebraic  equations  defining  the  a 

ij 

OL^  and  y. 

For  later  use  we  now  introduce  the  following  dimensionless  measures 

of  A , B und  a : 
mn  m n 


. _ 4^  , 2,2 

A = IT  Da  / pa  b 
mn  mn 


B H 2B  /pab 
m m 


A = 2a  /pab 
n n 


>(B22) 
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simplification  of  simultaneous  equations.-  From  (B18) 


Y - -pab/2 

With  this  result  substituted  Into  (B19),  (B20)  and  (B21) , the 
latter  become,  after  minor  algebraic  manipulation, 

Aio  - 0^B^-l)/i^  (1-  2,  4,  ....  M) 

- (Aj-l)/(j^B^')  (J  - 2,  4,  ...,  N) 

9 9 99  9 /1“2,4, 

A^j  - 2(Aj  + B^)3  (1  + - 2,  4 


These  can  be  used  to  eliminate  the  a In  equations  (B5)  and 

mn 

(B6) . Those  equations  can  then  be  written  as  follows: 


A « [n"^3“^  -I  2B  e^(m^  + n^3^)"^]^f  (n) 
" m=2,4...  “ ^ 


(n  - 2,  4 N) 


B f (m)  - m'^B^  - f 2A  B^(m^  + 
n-^  4 " 

11  • H j • • 

(m  - 2 , 4 , . . . , M) 


where 

f (n)  = n“^B"^  + I 2B^(m^  + n^B^)"^ 
m=*2,4, . • 


f (m)  = m ^B^  + I 2B^(ra^  4-  n^B^)”^ 

n-2,4, . . 


If  (B27)  is  now  used  to  eliminate  the  A in  (.B28) , the  latter 

n 

equations  become 
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(B23) 

(B24) 

(B25) 

(B26) 


(B27) 

(B28) 


(B29) 


M 

y B c r 


(m  - 2 , 4 , . . . , M) 


(B30) 


where 


N 


V 2 , 2„2.-2,  2 ^ 2„2.-2r-  , m 

- ~ ■ I 43  (m  + n 3 ) (q  + n 3 ) If,(n)] 

2 qm  ^.2.4,..  ^ 


-1 
(B31) 


c H m ^3^  - \ 2n  ^[f,(n)] 

m o / ■!■ 

e" 2 1 4 1 • • 


N 


-4,^2  . ^2„2,-2, 


-1 


On  the  basis  of  these  developments  we  now  have  the  following  cal- 
culation procedure:  Solve  equations  (B30)  (whose  number  depends  only 

on  M)  simultaneously  for  the  B^,  Wi ’-.h  these  known,  equations  (B27)  will 
furnish  the  values  of  the  A^.  Then  equations  (B24)  to  (B26)  will  yield 
all  the  except  A^^.  The  latter  is  defined  by  equation  (B4)  as 


00 


1 N M 

“ I K + I 


9 / *0n  ’ h,  , ^mO^ 

n“2,4, . . m«2,4, . . 


(B32) 


Deflections  and  edge  bending  moments,-  Once  the  A^j  are  known, 
the  deflections  and  edge  bending  moments  are  readily  computed.  Equation 
(Bl)  gives  the  former  as 


w(x,y)  = 


M 


N 


. I I A cos  cos  ^ (B33) 

.V  m-0.2,..  n-0,2...  ^ ^ 


The  bending  moments  normal  t<5  the  edge  x 0 are  also  obtained  with 
the  aid  of  equation  (Bl)  as  follows: 

N M 


-D(a 


2 -2 


a w/8x  =>  pb  IT  y y ^ cos(niry/b) 

n"0 , 2 , . . m*»  2 , 4 , . . 

(B34) 
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Those  along  the  edge  y ■ 0 are  obtained  In  a similar  way: 


(M  ) „ “ -D(3^w/3y^)  « “ pa^it  ^ \ T A n^cos(mTrx/a) 

(B35) 

Vertical  reactions  along  the  edges.-  The  unknown  vertical  reactions 
(force  per  unit  of  length)  along  the  edges  x ■ 0 and  y 0 will  be 
sought  In  the  following  Fourier  series  form: 

N 

(V  ) „ “ na  y S coa(mry/b) 

n-0,2...  ^ 

(B36) 

M 

^'^v^v-0  " 5!  T cos(m7Tx/a) 

where  the  coefficients  S and  T are  as  yet  unknown  except  for  the  fact 

n m 

that  overall  equilibrium  of  the  plate  requires  that 
pab  - 2pa(Syb)  + 2pb(TQa) 

from  which  it  follows  that 

Sq  + Tq  - 1/2  (B37) 

In  order  to  determine  the  and  T^  we  shall  postulate  an  auxilary 
problem,  that  of  a uniformly  loaded  plate,  with  edges  constrained  to 
have  zero  normal  slope,  with  prescribed  vertical  loadings,  in  the  form 
of  (B36) , along  the  edges  x - 0,  y - 0,  and  with  the  loadings  along  the 
other  two  edges  dictated  by  symmetry.  We  shall  determine  the  deflection 
shape  w(x»y)  for  this  p].ate  by  means  of  the  energy  method,  and  then  ask 
what  the  prescribed  vertical  loadings  along  the  edges  have  to  be  in 
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order  to  reduce  the  edge  deflections  to  zero,  l.e.,  to  bring  the  plate  of 
the  auxiliary  problem  into  the  same  state  of  deformation  as  the  plate 
of  the  original  problem. 

To  solve  the  auxiliary  problem,  we  shall  again  assume  the  deflections 
In  the  form  of  the  series  (Bl) , which  again  leads  to  equation  (B7)  for 
the  strain  energy  V^.  The  potential  energy  of  the  applied  loads  must 
now  Include  the  potential  energy  of  the  prescribed  vertical  loadings 
along  the  edges.  Thus,  in  addition  to  (B8),  we  have  the  following 
contribution  to  this  energy*. 

a b 

■ 2/  <V  + 2/  (V  ) w(0,y)dy  (B38) 


Substitution  of  the  aeries  expansions  (B36)  and  (Bl)  converts  this  tc 


M N T S 

V4  - 2pab  I I + 2^- 

m*0 , 2 , . . n*0 , 2 , . . mO  nO 


(B39) 


The  total  potential  energy  (TFE)  can  now  be  written  as 


TPE  - V,  + V-  + V, 
13  4 


(B40) 


where  the  three  terms  on  the  right  are  defined  by  equations  (B7), 

(B8)  and  (B39),  respectively. 

To  find  the  deflection  shape  under  the  given  loading,  we  now  make 
the  TPE  stationary  with  respect  to  che  a^j  via  the  equations. 


3(TPE)/3aQQ  - 0 


3(TPE)/3a^Q-  0 (i  - 2,  4*  . . . , M) 
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^(841) 

(Cont'd  on  next  pg.) 


3(TPE)/3a^,  - 0 (j  - 2,  4,  ....  N) 

Oj 

3(TPE)/Da^,  - 0 (1  - 2,  4,  ....  M;  j - 2,  4,  N) 

With  the  differentiations  carried  out,  these  become 

-pab  + 2pabTQ  + 2pabSQ  ■ 0 (B42) 

4 4 ’ 

* 4a,„  + pabT,  + 2pabS_  - 0 (i  - 2,  4 M)  (B43) 

Bab  2 iO  i 0 


Bab 


• 4a^^  + 2pabTQ  + pabS 


Oj 


'j 


(j  - 2,  4, 


N) 


(B44) 


Bab 


+ j^B)^*2a^j  + pabT^  + pabS^ 


« i “ 2,  4 M 

° ( ) 

- 2,  4 

(B45) 


At  this  point  the  similarity  between  this  group  of  equations 
and  equations  (BIB)  to  (B21)  should  be  noted. 

We  now  postulate  that  the  loading  functions 
arbitrary,  but  are  whatever  is  necessary  to  give  zero  deflection  along 
the  edges.  In  that  case  the  plate  deflections  of  the  auxiliary  problem 

must  be  identical  with  those  of  the  original  problem.  That  is,  the  a^^ 
in  equations  (B42)  to  (B45)  will  be  the  same  as  the  a^j  in  equations 
(BIB)  to  (B21).  It  then  follows  that 
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2pab(TQ+SQ>  - -2y 


pab(T^+2SQ>  - -(y+3^) 

(B46) 

pab(2TQ+Sj)  ■ “(Y+aj) 
pab(T^+Sj)  --(a^+e^) 

Substituting  y « -pab/2  (from  eq.  (B23)),  dividing  through  by  pab/2, 
and  making  use  of  tho  definitions  (B22) , we  can  convert  these  equations 
to  the  following  form: 


4(Tq+So)  - 2 

(B47) 

2(T^+2Sq)  - 1 - B^ 

(B48) 

2(2Tq+Sj)  - 1 - Aj 

(B49) 

2(T^+Sj)  - -(Aj+B^) 

(B50) 

where  1 ■ 2,  4,  ....  M and  J « 2,  4,  N.  The  first  of  these 

equations  confirms  the  equilibrium  condition,  (B37);  and  the  fourth 
equation  is  implied  by  the  first  three. 

We  thus  have  arrived  at  the  following  relationships  between  the 
Fourier  coefficients  in  equations  (B36)  and  the  dimensionless  Lagrangian 
multipliers  and  B^  defined  by  equations  (B30)  and  (B27)*. 


0 ■ 

(B51) 

(i  - 2.  4. 

....  M) 

(B52) 

(j  - 2,  4, 

....  N) 

(B53) 
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We  see  that  the  and  have  been  determined  to  within  a single 
undetermined  constant,  say  Sq.  The  removal  of  this  indeterminacy  is 
discussed  in  the  next  section. 

Removal  of  indeterminacy.  - We  note  first  that  the  choice  of  a value 
for  Sq  can  in  no  way  affect  the  deflections,  inasmuch  as  they  are  com- 
pletely determined  by  the  and  through  equations  (B24)  to  (B26) , 
(B32)  and  (B33).  On  the  othct  hand,  the  value  chosen  for  will  affect 
the  and  (see  eqs . (B51)  to  (B53))  and,  therefore,  the  reactions 


y y 

The  question  naturally  arises,  how  can  variations  in  Sq  alter  the 
vertical  reactions  along  the  edges  twhlch  can  iust  as  well  be  regarded 
as  loads)  without  at  the  same  time  altering  the  deflections?  A plausible 
conjecture  one  can  make  in  order  to  answer  this  question  is  that  equations 
(B36)  represent  not  only  distributed  reactions  but  also  a set  of  concen- 
trated self-equilibrating  forces  infinitesimally  close  to  the  corners,  us 
shown  in  Figure  15,  and  that  by  varying  we  are  varying  the  magnitude 
of  these  forces.  This  conjecture  is  plausible  because  concentrated  forces 
such  as  shown  in  Figure  15  will  alter  the  mean  values  of  (V  ) „ and 

and,  therefore,  and  Tq,  without  producing  any  deflections  and 
without  changing  the  sum  + T^. 

Accepting  this  conjecture,  it  now  becomes  clear  that  should  be  so 

chosen  as  to  cause  that  component  of  the  reactions  represented  in  Figure 

15  to  vanish.  The  reason  is  that  these  concentrated  forces  correspond  to 

Dirac  delta  function  components  in  (V  ) „ and  (V  ) and  if  the  series 

. X x-0  y y“0’ 

(B36)  have  to  represent  such  components  they  will,  not  converge  in  the 
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uaual  sense,  because  their  coefficients  S and  T will  not  approach  zero 
as  N,  M,  n and  m approach  infinity.  Thus,  if  the  reaction  component  shown 
In  Figure  15  is  not  "swept  out"  of  equations  (B36)  those  equations  will  be 
unusable  for  computing  anywhere  along  the  edge  x ■ 0 or  anywhere 
along  the  edge  y =t  0 (not  only  at  the  corners). 

It  has  already  been  mentioned  that  the  presence  of  any  P other  than 
zero  in  the  reactions  along  the  edges  will  manifest  itself  through  and 
not  approaching  zero  as  N,  M,  n and  m approach  infinity.  Thus,  either 
of  the  following  conditions  should  suffice  to  eliminate  any  P-forces  from 
equations  (B36) : 


Llm 

Llm 

M-x» 


"m-  0 


(B54) 


Within  the  framework  of  a solution  with  a finite  number  of  terms,  and 
T|^  themselves  are  the  best  available  estimates  of  the  left-hand  sides. 
Thus,  equations  (BSA)  will  be  replaced  by 


(B55) 


and.  In  view  of  equations  (B52)  and  (B53),  these  lead  to  the  following 
equations  for  T^  and 


Sq  - ~ Bm>  (®57) 


The  following  two  calculation  procedures  now  suggest  themselves: 

(a)  Use  (B57)  to  compute  S^,  then  (B51)  to  get  T^,  and  accept  some  degree 


-47- 


of  non-satisfaction  of  (B56).  (b)  Use  both  (B56)  and  (B57)  to  evaluate 

Tq  and  Sq  and  accept  some  degree  of  non-satisfaction  of  (B51) . In  the 
second  procedure  the  closeness  of  Tq  + Sq  to  1/2  can  be  taken  as  an  over- 
all measure  of  the  convergence.  The  second  procedure  Is  the  one  that  was 
used  In  the  computations  leading  to  the  K - «>  (arctan  K ■ ’n/l)  curve  of 
Figure  6. 

An  alternate  hypothesis  was  also  Investigated  as  a basis  for  removing 
the  Indeterminacy  In  Sq.  This  Is  the  hypothesis  that  at  a corner  the 
vertical  reaction  Intensity  Is  continuous  around  the  corner.  That  Is, 


*''^x^x"0,y-0  " ^'^y^y-0,x-0 


(B58) 


or,  in  view  of  (B36) , 
N 


M 


a I S - b ^ T 
n m 

n^O.Z, . . m-0,2, . . 


(B59) 


Eliminating  and  through  (B52)  and  (B53)  , then  replacing  Tq  by  y - Sq, 
and  then  solving  for  Sq,  one  obtains 

N M 

2+BN+M+2ey  k - 2 ) B 

^ n ^ m 

g _ n»2,4, . . m«2,4,. . 


4(1+N)  (1+6) 


(B60) 


as  the  equation  defining  Sq  on  the  basis  of  this  second  hypothesis.  The 
alternate  procedure  just  described  for  removing  the  Indeterminacy  in  Sq 
leads  to  essentially  the  same  numerical  results  as  the  first  procedure  . 
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APPENDIX  C 


REMARKS  ON  THE  CORNER  STRESS 


The  analysis  of  appendix  A,  which  produced  the  values  of  R given  in 
Figure  7,  is  based  on  the  assumption  that  the  rotational  spring  stiffness 
k (see  '’0.  (1))  furnished  by  a waxl  to  the  adjacent  odge  of  the  lid  is 
constant  all  along  the  length  of  that  edge.  Since  each  wall  is  essent- 
ially a wide  beam  this  is  a plausible  assumption  everywhere  except  near 
the  corners.  Where  two  walls  meet  at  a corner  their  continuity  with  each 
other  tends  to  increase  the  rotational  stiffness  of  each;  in  fact,  this 
stiffness  approaches  infinity  right  at  the  corner.  This  discrepancy  between 
the  actual  iiature  of  the  rotational  restraint  and  that  assumed  in  the 
analysis  is  purely  a local  one.  As  long  as  the  height  h of  the  package  is 
small  compared  to  the  dimensions  a and  b there  should  be  very  little  error 
in  any  computed  quantities  except  for  the  corner  force  R.  R is  proportional 
to  the  race  of  twist  of  the  lid  at  the  corner  and  that,  in  turn.  Is  strongly 
depeiiJent  on  the  local  constraint  conditions  at  the  corner.  Consequently, 
it  is  highly  probable  that  the  R values  furnished  by  Figure  7 are  not  valid 
(probably  too  high)  when  applied  to  the  lid  of  a flat-pack  under  external 
pressure. 

Other  considerations  lead  to  the  conclusion  that  the  corner  reactions 

for  the  flat-pack  lid  might  even  be  zero:  If  one  assumes  perfect  joint 

rigidity  at  the  corner,  not  only  between  the  two  walls  but  also  between 

each  wall  and  the  lid,  it  then  follows,  by  a not  too  difficult  argument, 

that  the  rate  of  twist  of  the  lid  must  vanish  at  the  corner.*  Consequently, 

*Referring  to  Figure  13,  consider  three  points  A,  B and  C infinitesimally 
close  to  the  corner  0.  Assume  perfect  joint  rigidity,  which  implies  that 
right  angles  are  maintained  at  points  1,  2 and  3 between  the  pair  of  dotted 
lines  meeting  at  each  of  these  points.  Now  suppose  there  is  some  twist  of 
the  lid  at  the  corner.  This  twist  will  cause  point  A to  move  vertically, 
say  downward.  The  mainte  ance  of  right  angles  at  points  1 and  2 will  then 
require  that  points  B and  t both  move  outward.  However,  this  would  destroy 
the  right  angle  which  must  exist  at  point  3 between  lines  C3  and  B3.  We 
must,  therefore,  conclude  that  point  A can  have  no  vertical  movement;  i.e., 
the  lid  can  have  no  twist  at  the  corner. 
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R,  which  is  proportional  to  the  rate  of  twist  at  the  corner,  must 
vanish. 

We  are  thus  faced  with  an  uncertainty  regarding  the  magnitude  or 

even  the  existence  of  the  tensile  stress  in  the  seal  at  the  corners  of 

the  package.  Fortunately,  for  most  packages  the  proportions  are  such 

that  as  given  by  equation  (12),  will  exceed  given  by 

equation  (19).  For  those  packages  will  not  govern  in  either  the 

design  or  the  screening  process,  and  the  question  as  to  whether  it  exists 

or  whether  equation  (12)  over-estimates  it  becomes  academic.  For  those 

rare  cases  in  which  S , as  given  by  equation  (19),  might  govern,  a 

corner 

conservative  approach  to  this  question  has  been  outlined  in  the  body  of 
this  report. 
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Figure 


(a) 


Figure 


1.-  Flat-pack  configuration. 
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(b)  Bottom  edge  subjected  to 
same  moment  as  upper  edge 


2.  Moment  rotation  relationships  for  walls. 
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Figure  3. 


arc tan  K 


- Relationship  between  K,  K 


1 


and  arctan  K. 
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.2 


arctan  K 

Figure  5.-  Bending  moment  intensity  at  the  middle  of  the  long  side  for 
a uniformly  loaded  rectangular  plate  with  edges  elastically 
restrained  against  rotation. 
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Figure  7.-  Magnitude  of  concentrated  corner  reaction  for  a uniformly 
loaded  rectangular  plate  with  edges  elastically  restrained 
against  rotation.  (Applies  for  - 0.3.  For  other  values 
of  V multiply  ordinate  by  (l-v^/0.7.) 
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Figure  8. 


arctan  K 

- Center  deflection  of  a uniformly  loaded  rectanguJLar  plate  with  edges 
•laatlcally  restrained  against  rotation  (small-d*f lection  theory). 
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(a)  Uue  to  M (b)  Due  to  V 

' max  max 

Figure  1.0.-  Streau  distributiona  across  the  width  of  the 
seal  at  the  middle  of  the  long  side. 
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Figure  11.-  Area  (shaded)  over  wliich  the  corner  reaction  R is 
assumed  to  be  distributed. 
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Area  (shaded)  over  which  corner  reaction  R is  assumed 
to  be  distributed  for  package  with  rounded  corners 
where  the  walls  meet. 


Points  near  the  Junction  of  two  walls  and  the  lid. 
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(a) 


(b) 


Figure  14.-  Two  possible  kinds  of  lid-wall-base  interaction 
for  a package  in  a centrifuge. 


Figure  15.- 


A loading  which  will  not  produce  deflections  but 

and  (V  ) 


will  alter  the  Fourier  series  for  (V^) 


x=0 
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